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We theoretically study the population inversion and negative dynamic conductivity in intrinsic 
graphene in the terahertz (THz) frequency range upon pulse photoexcitation with near- /mid- infrared 
wavelength. The threshold pulse energy required for the population inversion and negative dynamic 
conductivity can be orders-of-magnitude lower when the pulse photon energy is lower, due to the 
inverse proportionality of the photoexcited carrier concentration to the pulse photon energy and 
to the weaker carrier heating. We also investigate the dependence of the dynamic conductivity on 
the momentum relaxation time. The negative dynamic conductivity takes place either in high- or 
low-quality graphene, where the Drude absorption by carriers in the THz frequency is weak. 
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I. INTRODUCTION 

Graphene has attracted much attention for wide vari- 
ety of device applications due to its exceptional electronic 
and optical properties. 1 3 Especially, THz devices such as 
lasers 4 9 and photo detect or s2&~— utilizing its high carrier 
mobility and gapless dispersion have been investigated 
(for a review of these devices, see Ref.IH). In Refs.HM 
we demonstrated that the population inversion can occur 
in optically pumped graphene at THz/far-infrared range 
of frequency and hence the lasing at such a range is pos- 
sible, utilizing the gapless linear energy spectrum and 
relatively high optical-phonon (OP) energy in graphene. 
Due to the energy spectrum e = ±v F hk, where v F c± 10 8 
cm/s is the Fermi velocity and k is the wavenumber, the 
Fermi energy e F in equilibrium in intrinsic graphene is 
equal to zero. Hence, the electron and hole distribution 
functions at the bottom of the conduction band and the 
top of the valence band have values / e (0) = /h(0) = 1/2. 
This implies that upon photoexcitation one might expect 
to make the values of the distribution functions at low 
energies greater than one half, i.e., f e (e) = fh(z) > 1/2, 
corresponding to the population inversion. Such a popu- 
lation inversion leads to the total dynamic conductivity in 
graphene at THz/far-infrared frequencies being negative, 
and lasing in graphene at such frequencies is possible. 

Recently, we measured the carrier relaxation and 
recombination dynamics in optically pumped epitax- 
ial graphene on silicon 14 and exfoliated graphene 15 us- 
ing THz time-domain spectroscopy based on an optical 
pump/THz&optical probe technique, and observed the 
amplification of THz radiation by stimulated emission 
from graphene under pulse excitation. To our knowl- 
edge, those are the first observation of the THz amplifi- 
cation using optically pumped graphene. Those results 
demonstrate the possibility of realizing THz lasers based 
on graphene. 

The carrier dynamics in optically pumped graphene 
strongly depends on the initial temperature of carriers 
and the intensity of the optical pumping. With suffi- 



ciently low carrier concentration, i.e., at low tempera- 
ture and upon weak pumping, photoexcited carriers are 
effectively accumulated near the Dirac point via cas- 
cade emission of OPs, and the achievement of the popu- 
lation inversion is expected to be efficient On the 
contrary, at room temperature or upon not so weak 
pumping where the carrier concentration is rather high, 
the carrier-carrier (CC) scattering plays a crucial role 
in the dynamics after the pulse excitation, due to fast 
quasi-equilibration of carriers. In fact, ultrafast optical 
pump-probe spectroscopy on graphene has indicated that 
the quasi-equilibration by the CC scattering takes place 
within the time scale of 10 fs^~— which is much faster 
than the process of single OP emission. In such a case, 
the pulse excitation makes carriers very hot initially, and 
the energy relaxation and recombination via OP emission 
follow. 

Previously, we studied the population inversion and 
the negative dynamic conductivity in intrinsic graphene 
at room temperature upon the pulse excitation with pho- 
ton energy 0.8 eV, corresponding to the wavelength 1.55 
/im, and showed that they can be achieved with the pulse 
energy fluence above the thresholdi2£ 

In this paper, we investigate the dependence of the 
threshold pulse fluence for the THz negative dynamic 
conductivity on the pulse photon energy using the model 
extended from that developed in Refs. @ and [2(J to take 
into account the effect of the Pauli blocking. We also 
study the dependence of the dynamic conductivity on 
the momentum relaxation time r. We consider the two 
limiting cases where ujt > 1, corresponding to high- 
quality graphene, and where ujt <C 1, corresponding to 
low-quality graphene. 



II. EQUATIONS OF THE MODEL 

Assuming the quasi-Fermi distribution of carriers due 
to the quasi-equilibration by the CC scattering, i.e., 
f e = {1 + exp[(£ — £f)/^c]} _1 , the time evolution of 



2 



the distribution is represented through the quasi-Fermi 
level £p(t) and the carrier temperature T c (t). Here, both 
electron and hole distributions can be expressed by a sin- 
gle "carrier distribution" because we consider intrinsic 
graphene; due to the symmetric dispersion for electrons 
and holes in the energy range under consideration, i.e., 
e < 1 eV, their distributions remain identical even upon 
the pulse excitation. The rate equations which deter- 
mine the quasi-Fermi level and carrier temperature can 

be obtained from the quasi-classical Boltzmann equa- 
tions^ 
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where n = n(eF,T c ) and £ — £(sf,T c ) are the concentra- 
tion and energy density of either type of carriers, which 
are found by integrating over k the distribution function 
multiplied by proper factors, the index i runs over types 
of OPs (long- wavelength T-OP and short- wavelength K- 

OP), and Tif^tra an d r f inter are ^ ne intraband and inter- 
band scattering rates for OPs ((+) for absorption and 
(— ) for emission). Those rates are on the order of sub- 
picoseconds for high-energy carriers. 

Note that the first term (intraband transition) in the 
curly brackets in Eq. ([2]) contains the factor fwji , whereas 
the second term (interband transition) contains e. These 
are energy that a carrier looses/acquires in their corre- 
sponding transitions. The energy conservation for the 
latter transition is implicitly satisfied because the rest of 
the OP energy huji — e is taken into account in the same 
equation for the other kind of carrier. Here, we assume 
the equilibrium OPs and neglect the nonequilibrium pop- 
ulation of OPs due to their emission via OP scattering 
of carriers. The effects of the latter have been discussed 
elsewhere^ 

Equations (pQ) and (|2]) together with initial conditions 
form a nonlinear system of equations for £ F (t) and T c (t), 
which can be solved numerically. Equations (pQ) and (|2]) 
are accompanied with the initial carrier concentration 
and energy density from which the initial quasi-Fermi 
level and carrier temperature can be found: 
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where no and £q are the intrinsic carrier concentration 
and energy density, and An and A£ are contributions of 



photogenerated carriers. Those are equal to the concen- 
tration and energy density generated from the pumping: 

iraAJ . p . <. A „ HQ A 
An~-^(l-2/ Ml/2 | t=0 ), A£~— An, (4) 

where a ~ 1/137, e is the dielectric constant surrounding 
graphene, Ml and AJ are the pulse photon energy and 
the pulse fluence. 

In Eq. the last factor 1 — 2/^/ 2 roughly takes 
into account the Pauli blocking of the photoabsorption, 
and the distribution function in it is approximated by 
the quasi-Fermi distribution after the photoabsorption 
with SF\t=o and T c \ t =o calculated self-consistently from 
Eqs. (|3]) and (jU). This approximation underestimates 
the effect of the Pauli blocking by the nonequilibrium 
photogenerated carriers when the pulse is so intense that 
the rate of the photogeneration is comparable to that of 
the quasi-equilibration by the carrier-carrier scattering. 
On the other hand, it overestimates the effect by time- 
dependent increase in the distribution function during 
the pulse excitation, somewhat compensating the under- 
estimation described above. In addition, we assume the 
pulse has sufficiently short width that the OP scattering 
during the absorption of the pulse is not effective, i.e., 
Ai<^e,W2),i» ra (ffi/2). 

To study not only the gain acquired by the popula- 
tion inversion but also the net gain taking into account 
the loss due to the Drude absorption of electromagnetic 
waves by carriers in graphene, we introduce the real part 
of dynamic conductivity 22 as it is related to the net gain: 
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where r is the momentum relaxation time. The first 
term in Eq. (j5j) corresponds to the interband contribu- 
tion which depends on the frequency only through the 
distribution function and which can be negative when 
the rate of stimulated emission exceeds the rate of ab- 
sorption, i.e., when the population inversion takes place. 
Since the distribution function is the quasi-Fermi distri- 
bution, the condition of the population inversion for the 
photon energy two is simply represented as Sp > hoo/2. 
On the other hand, the second term corresponds to the 
intraband contribution (Drude conductivity) which is al- 
ways positive in the system under consideration. 

When the real part of the dynamic conductivity be- 
comes negative at some frequency, the electromagnetic 
wave at that frequency passing through optically pumped 
graphene is amplified and it serves as a gain medium. 
Note that from Eq. ([5]) the largest achievable negative 
value of the conductivity is — e 2 /4ft, corresponding to the 
amplification of (2.3/ y/e) % of the incoming electromag- 
netic wave. 
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FIG. 1. The time evolution of quasi- Fermi level and carrier 
temperature with different pulse fluence and pulse photon en- 
ergy [(a) HQ = 0.8 eV, (b) 0.4 eV, (c) 0.2 eV, (d) 0.1 eV]. 
Room temperature (300 K) is indicated by thin solid lines. 



III. RESULTS AND DISCUSSION 

A. Time Evolution of Quasi-Fermi Level and 
Carrier Temperature 

Using the model above, we investigate the population 
inversion and the time-dependent dynamic conductivity 
in intrinsic graphene at room temperature after pulse 
excitation. We varied the pulse photon energy HQ in 
the range between 0.1 — 0.8 eV corresponding to infrared 
wavelength of 12.4 — 1.55 fim. Besides, we set e = 5.5 
corresponding to the effective dielectric constant of the 
interface between the air and SiC. Here, we have used the 
pulse fluence rather than the pulse intensity as a measure 
of pulse strength because the former does not depends 
explicitly on the pulse photon energy, whereas the latter 
does through the pulse width. 

Figure [T] shows the time evolution of the quasi- Fermi 
level and carrier temperature after the pulse excitation 
with different pulse fluence and pulse photon energy. It 
is clearly demonstrated in Fig. [I] that the population in- 
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FIG. 2. The dependence of the threshold fluence of the pop- 
ulation inversion on the pulse photon energy for different fre- 
quencies of THz photons. The inset shows the same depen- 
dence but with linear scale. 



version at THz frequencies, say, up to about 10 THz cor- 
responding to fiuj/2 ~ 20 meV is achieved with some 
threshold pulse fluence, and that it lasts on the order of 
10 ps. 

The mechanism of the population inversion in the THz 
range is explained as follows. Right after the pulse ex- 
citation, the carrier temperature becomes very hot (ex- 
cept for the excitation with low photon energy) and the 
quasi-Fermi level can even be negative due to the quasi- 
equilibration of carriers with very large heat brought by 
the photogenerated carriers. After a few ps, carriers with 
high energy are relaxed via intraband OP emission and 
accumulated in the low energy region, which is illustrated 
by the increase in the quasi-Fermi level. As discussed in 
Ref. [20|, the population inversion lasts on the order of 10 
ps due to the imbalance between the time scales of intra- 
and inter-band OP emission. 

Comparing Figs. [TJa)-(d), one can see that the thresh- 
old fluence of the population inversion becomes lower as 
the pulse photon energy becomes lower. This is also seen 
in Fig. [2] that illustrates the dependence of the threshold 
fluence of the population inversion on the pulse photon 
energy for different frequencies of THz photons. There 
are several factors for this lowering. First, the lower flu- 
ence is needed for the pulse with lower photon energy 
to have the same amount of carrier concentration [more 
precisely, A J oc M^, as seen in Eq. (H])]. Second, the heat- 
ing of carriers becomes less and the carrier temperature 
becomes lower [even the cooling can occur as shown in 
Fig. QJd); see below for more details], resulting in their 
denser accumulation in the low energy region. Finally, 
the interband OP emission, i.e., the recombination is sup- 
pressed due to the lower carrier temperature. All these 
result in the nonlinear lowering of the threshold as the 
pulse photon energy decreases; Fig. [2] exhibits almost lin- 
ear dependence of the threshold in the region of large 
photon energy, caused by the first factor, and the sharp 
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FIG. 3. The time-evolution of the dynamic conductivity in the THz range normalized to e 2 /Ah and cut off above zero. The 
columns from left to right correspond to cases with Ml = 0.8, 0.4, 0.2, and 0.1 eV and AJ = 6 x 10 3 , 1 x 10 3 , and 1 x 10 2 
nJ/cm 2 , respectively, whereas the rows from top to bottom correspond to cases with r = 10 ps, 1 ps, and 1 fs, respectively. 



drop in the region of photon energy below 0.25 eV, caused 
by the second and third factors. 

The effect of the Pauli blocking depends strongly on 
the pulse photon energy. For Ml = 0.8 eV, the changes in 
the quasi-Fermi level and carrier temperature were neg- 
ligible for any pulse fluence in consideration when the 
Pauli blocking was taken into account. On the contrary, 
for HQ = 0.1 eV, the threshold fluence changed 1.3 — 1.8 
times larger depending on the frequency, and the quasi- 
Fermi level achievable became smaller. 

As has been discussed in Refs. [B| and [H, the cooling 
in the case of low pulse photon energy occurs because 
the initial photoexcitation can give the lower energy to 
carriers than in thermal equilibrium. In fact, the en- 
ergy of photogenerated carriers, HQ/2, is lower than the 
average energy of carriers in intrinsic graphene at room 
temperature, S/n « 0.0567 eV, when HQ < 0.13 eV. 
Besides, the cooling after the heating by the photoexci- 
tation (albeit small) can also occur as seen in Figs. GJa)- 
(c). This is associated with the discrete nature of OP 
emission/absorption which cannot exactly bring the car- 
rier temperature to the lattice temperature. The latter 
is to be accomplished via the acoustic phonon scatter- 
ing and/or radiative generation/recombination with sub- 
stantially longer time scale. 



B. Time-Evolution of Dynamic Conductivity in 
Cases of High- and Low-Quality Graphene 

The dynamic conductivity in Eq. (j5j) depends primarily 
on the carrier distribution, i.e., the carrier concentration, 
as well as the momentum relaxation time r. Since the 
photoexcitation not only leads to the population inver- 
sion but also to the increase in the carrier concentration, 
the negative dynamic conductivity arises if the increase 
in the absolute value of the interband conductivity over- 
comes the increase in the Drude absorption in Eq. (|5]). 
For the realization of low threshold of negative dynamic 
conductivity in the THz range, it is necessary to suppress 
the Drude absorption. This can be done by using either 
good-quality graphene such that out > 1 or high-quality 
graphene such that ujt <C 1. As can be seen in Eq. (|5]), 
in the former case, the high-frequency "cut-off" of the 
Drude conductivity occurs. In the latter case, the Drude 
conductivity ceases in the entire frequency range. 

Figure [3] illustrates the time-evolution of normalized 
dynamic conductivity Re cr^ / (e 2 / 4H) in the THz range 
with different pulse photon energy and momentum re- 
laxation time. The pulse fluence for each pulse photon 
energy was chosen so that sufficiently large negative dy- 
namic conductivity is obtained. It demonstrates that the 
negative dynamic conductivity can be achieved in the 
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FIG. 4. The maximum negative value of the dynamic conductivity in its time duration as a function of frequency and pulse 
fluence. The columns from left to right correspond to cases with HQ — 0.8, 0.4, 0.2, and 0.1 eV, respectively, whereas the rows 
from top to bottom correspond to cases with r = 10 ps, 1 ps, and 1 fs, respectively. 



THz range with sufficiently high pulse fluence, and that 
the same level of negative dynamic conductivity can be 
achieved with much lower pulse fluence for lower pulse 
photon energy. It also shows that its magnitude, fre- 
quency range, and time duration depend strongly on the 
momentum relaxation time. 

For r = 10 ps, the negative dynamic conductivity oc- 
curs down to 1 THz and extends to more than 10 THz, 
and it duration is 5 — 30 ps depending on the frequency. 
The duration is longer for lower frequency as is naturally 
understood from the relaxation of the quasi-Fermi level 
shown in Fig. [TJ For r = 1 ps, the negative dynamic 
conductivity occurs between 4 and 10 THz with smaller 
magnitude and much shorter duration. This indicates 
that the high quality of graphene having longer momen- 
tum relaxation time is preferable. As can be seen in 
Fig. [3j the minimum frequency for the negative dynamic 
conductivity is primarily determined by the momentum 
relaxation time, although it slightly shifts to the lower 
side in case of the lower pulse photon energy. 

On the contrary, for r = 1 fs, the negative dynamic 
conductivity occurs in the low frequency region below 
6 THz. Note, however, that the meaning of the "time- 
dependent dynamic conductivity" is valid only if its vari- 
ation is slower than the frequency under consideration, 
approximately around or above a few hundred GHz in our 
case. This means we can have the negative dynamic con- 
ductivity down to a few hundred GHz. The upper limit 



of the frequency range having the negative dynamic con- 
ductivity in this case is smaller than in the case where 
lot <C 1. This is because the Drude conductivity in this 
case is almost constant in the THz range, whereas it van- 
ishes due to the high-frequency "cut-off" in the latter 
case. 

Our calculation showed that even for HQ = 0.1 eV 
the threshold momentum relaxation time of the nega- 
tive dynamic conductivity is r = 5 fs, and above it 
the dynamic conductivity does not occur in the THz 
range, no matter how large the pulse fluence is. This 
threshold roughly corresponds to the value at which the 
Drude term in Eq. ([5]) at zero-frequency limit exceeds the 
largest negative value of the interband contribution, i.e., 
T~(7r/81og2)(fi/* B T z ). 

Figure [H shows the maximum negative value of the dy- 
namic conductivity in its time duration as a function of 
frequency and pulse fluence with different pulse photon 
energy and different momentum relaxation time. It illus- 
trates that the threshold pulse fluence as well as the pulse 
fluence needed to obtain the same value of the negative 
dynamic conductivity are decreased in several orders of 
magnitude for any frequency as the pulse photon energy 
is decreased. Comparing the cases with r = 10 and 1 
ps, it is seen that the decrease in the momentum relax- 
ation time causes rather large increase in the minimum 
frequency for the negative dynamic conductivity and the 
increase in the threshold pulse fluence for frequency near 
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the minimum, although the latter for frequency suffi- 
ciently above the minimum is insignificant due to the 
high-frequency "cut-off" of the Drude conductivity. As 
mentioned above already, the minimum frequency for the 
negative dynamic conductivity shifts to the lower side for 
the lower pulse photon energy. In fact, the frequency can 
be down to about 700 GHz with Ml = 0.1 eV and r = 10 
ps, whereas it is around 1.1 THz with Ml = 0.8 eV and 
r = 10 ps, demonstrating the superiority of the pumping 
by lower photon energy. For r = 1 fs, the threshold pulse 
fluence for each photon energy is on the same order in 
comparison with that for r = 10 ps. 

Also, it should be noted that the maximum saturates 
below the "fundamental limit" (i.e., — e 2 /4h) for very 
high pulse fluence. For example, for r = 10 ps and HQ = 
0.1 eV, it saturates with the value Re cr^ / (e 2 / AH) ~ —0.4 
at / = 1 THz and « -0.8 at / = 2 THz. The rea- 
son for this is that, for very high pulse fluence, the in- 
crease in the Drude absorption due to the increase in the 
carrier concentration exceeds the interband contribution 
even though very high quasi-Fermi level is obtained, and 
the maximum negative value of the dynamic conductiv- 
ity corresponding to the saturation value is reached at 
later time after the relaxation of the carrier concentra- 
tion. The Pauli blocking also takes part in the saturation 
in cases of lower photon energies. The saturation values 
for r = 1 fs are smaller than those for r = 10 ps, espe- 
cially, about 3 times smaller with Ml = 0.4 and 0.8 eV. 
Again, this is due to the nonzero, almost-constant value 
of the Drude conductivity in the THz range in the limit 

OUT <C 1. 



IV. CONCLUSIONS 



We studied the population inversion and negative dy- 
namic conductivity in intrinsic graphene upon pulse pho- 
toexcitation with the pulse photon energy in 0.1 — 0.8 eV 
and their dependences on the pulse fluence. We showed 
that their threshold pulse fluence can be lowered by two 
orders of magnitude with the lower pulse photon energy. 
This is attributed to the inverse proportionality of the 
photogenerated carrier concentration to the photon en- 
ergy, as well as to the weaker carrier heating which results 
in the denser accumulation of carriers in the low energy 
region and in the suppression of recombination via in- 
terband OP emission. We also showed that the negative 
dynamic conductivity takes place in the long relaxation 
time limit ojt > 1, where the Drude absorption in the 
THz range becomes small by the high-frequency "cut- 
off", and in the very short relaxation time limit ojt <C 1, 
where it ceases in the entire frequency range. We demon- 
strated that the pulse excitation by lower pulse photon 
has both lower threshold pulse fluence and the lower min- 
imum frequency for the negative dynamic conductivity 
down to about 700 GHz with HQ = 0.1 eV and r = 10 
ps. 
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